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ABSTRACT 

We describe a 22-year survey for variable and transient radio sources, performed with 
archival images taken with the Molonglo Observatory Synthesis Telescope (MOST). 
This survey covers 2775 deg 2 of the sky south of 5 < —30° at an observing frequency 
of 843 MHz, an angular resolution of 45 x 45csc|<5| arcsec 2 and a sensitivity of 5c ^ 
14mJy beam -1 . We describe a technique to compensate for image gain error, along 
with statistical techniques to check and classify variability in a population of light 
curves, with applicability to any image-based radio variability survey. Among radio 
light curves for almost 30000 sources, we present 53 highly variable sources and 15 
transient sources. Only 3 of the transient sources, and none of the variable sources 
have been previously identified as transient or variable. Many of our variable sources 
are suspected scintillating Active Galactic Nuclei. We have identified three variable 
sources and one transient source that are likely to be associated with star forming 
galaxies at z ~ 0.05, but whose implied luminosity is higher than the most luminous 
known radio supernova (SN1979C) by an order of magnitude. We also find a class of 
variable and transient source with no optical counterparts. 

Key words: radio continuum: general, stars: pulsars:general, stars: flare, techniques: 
image processing, catalogues, supernovae: general 



1 INTRODUCTION 

Many of the most extreme events in the universe mani- 
fest themselves through variability. For example, explosive 
events such as Gamma Ray Bursts (GRBs) are thought 
to originate from stellar collapses which can momentar- 
ily o ut-shine the entire universe (|MacFadven fc Wooslevl 
1999); the magnetosphere of the Crab pulsar can pro- 
duce nanosecond- duration pulses with a brightness temper- 
ature of 10 38 K jHankins et all [2003); thermonuclear ex- 
plosions up to a month long have been observed in X- 
ray binaries l|Sokoloski et aljfeoOrj , iBode et al.l I2006T ): and 
large variations in accretion rate have been observed in ac- 
creting black holes, ranging in size from Galactic micro- 
quas ars dFender et a l. 1999), to supermassive black holes 
(e.g. iGezari et al.ll2006l . iBlazeiowski et al1l2005l 1. 

In addition to transient emission indicating extreme 
events, many aspects of the interstellar and intergalactic 
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media, and of the universe as a whole, can be revealed 
through measurements in the time domain. For example, 
the p resence of ma ssive compac t halo objects |Alcock et al.l 
Il997h and planets (|Gouldll2005h can be inferred from grav- 
itational microlensing of starlight around massive objects; 
the expansion history of the universe ca n be measured us- 
ing light curves of Ty pe la supernovae |Riess et alj|l99§ , 
iPerlmutter et af]|l999h ; and interstellar scintillation (ISS) 
of compact radio sources can be used to infer the distance, 
velocity and turbu lence of otherwise in visible ionised clouds 
in the Milky Way |Bignall et al.ll2003h . 



Low-energy radio photons can be easily produced by 
both thermal and non-thermal processes, and in contrast 
to optical and X-ray wavelengths, the universe is essentially 
transparent to these photons. As a vast range of objects 
produce radio emission, and these signals are detectable over 
much of the observable universe, radio variability surveys are 
an effective method of discovering and characterising such 
sources and propagation phenomena. Furthermore, impor- 
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tant propagation effects such as interstellar scintillation are 
only observed in the radio. 

While radio observations are often used for follow-up 
of detections at other wavelengths, blind radio surveys have 
also been a fruitful method fo r discovering new objects and 
phenomena in their own right. iGreeorv fc Tavlorl (|l986h dis- 
covered a number of flaring radio stars in a survey of the 
Galactic plane, while I Bower et all (|2007l l with 944 images 
of a single field, discovered ten faint tra nsient source s, eight 
of which had no optical counterparts. ICarilli et al.l (120031 ) 
also discovered a number o f highly variable sourc es with ob- 
servations of a single field. iLevinson et al.l ([2002) covered a 
large fr action of the sky w ith two epoch s by comparing the 
NVSS (jCondon et al.lll998l ) and FIRST ijBecker et al.lll995l ) 
images, and was able to s et con straints on the GRB beam- 
ing fraction. ICroft et al.l |2010l ) set an upper limit on the 
rate of transients with flux densities greater than 40 mjy 
of 0.004 deg~ 2 by comparing a 690 square degree image ob- 
tained with the Allan Telescope Array (ATA) at 1.4 GHz 
with NVSS. 

A n umber of transient sourc e detections remain unex- 
plained. iMatsumura et all (|2009h have discovered a num- 
ber of > Uy transient sources with ~ 1 day time-s cales 
at a range of Galactic latitudes. lLorimer et all (|2007l ) dis- 
covered a 30 Jy transient source at 1.4 GHz with a 5 ms 
time-scale, but there is ev id ence this may be atmospheric 
l|Burke-Spolaor et alj l20ich . iLenc et al l (|2008h found a ~ 
100 mjy transient at 330 MHz during deep Ver y Long Base- 
line In terferometer (VLBI) observations, and [Becker et al.l 
(|2010h found 39 variable sources in the Galactic Plane with 
3 epochs of 5 GHz Very Large Array (VLA) observations, 
most of which had no known counterparts at other wave- 
lengths. 

In this paper we present a survey for transient and vari- 
able sources at 843 MHz with characteristic time-scales from 
days to years. We aim to characterise the variable and tran- 
sient radio sky at this frequency, enumerate the most ex- 
treme variable sources, find transient sources and develop 
techniques suitable for upcoming radio surveys. 

For our analysis we are using the Molonglo Obser- 
vatory Synthesis Telescope (MOST). The MOST and its 
predecessor the Molonglo Cross, have been at the fore- 
front of research in this area o ver se veral decades. Using 
the Molonglo Cross, iHunsteadl (|l972l ) was the first to ob- 
serve low frequency variability at 40 8 MHz. In an archival 
survey of calibrator measurements, iGaensler fc Hunsteadl 
(2000) found that one-third of the bright point sources 
at 843 MHz were variable, and a weak Galactic latitude 
dependence was demonstrated, indicating that interstellar 
scintillation was at least partly responsible. The MOST 
has been used to d iscover several hu ndred pulsars in a 
number of surveys l|Manchesterl Il985l ). and a survey for 
short- duration transients has been performed by I Amy et al.l 
|l989l) but with a null result. As a follow-up instrument, 
the MOST was the first telescope to dete ct prompt radio 
emission from SN1987A (|Turtle et al.lll987h . and it has been 
used to monitor a num ber of Galactic accreting systems 



Hannikaincn et al. 1998), a brightening supernova remnant 



Murphy et ail [2008) and a magnetar flare (|Gaensler et alj 



120051 ). 

The MOST archive is a unique resource which is able 
to address the limitations of sky coverage, sensitivity and 



cadence that have accompanied other blind surveys. In ad- 
dition, the experience gained from analysing wide field- 
of-view images provides an opportunity to develop tech- 
niques suitable for upcoming wide-field transient and vari- 
ability surveys such as the Variables And Slow Transients 
(VAST) survey for th e Australian Square Kilomet re Array 
Pathfinder (ASKAP) (IChatteriee fc Murphvll2010l) and the 
LOFAR Transients Key Project (jFender et al.ll2008l ). 

In §2 we describe the MOST and its image archive. In 
§3 we describe our method of extracting light curves from 
this archive and in §4 we present the results of applying our 
method, including quality checks and selected sources. In § 5 
we discuss our results, and in §6 we discuss some noteworthy 
sources. We draw our conclusions in §7. 



2 OBSERVATIONS WITH THE MOLONGLO 
OBSERVATORY SYNTHESIS TELESCOPE 

The Molonglo Observatory Synthesis Telescope (MOST) is 
located near Canberra, Australia and was constructed by 
modification of the East- West arm of the former One-Mile 
Mills Cross telescope. The MOST is an east-west synthesis 
array comprising two cylindrical paraboloids each of dimen- 
sion 778 m x 12 m separated by 15 m. Radio waves are 
received by a line feed system of 7744 circular dipoles. The 
telescope is steered in the North-South axis by mechanical 
rotation of the paraboloids about their long axis, and in the 
East- West axis by phasing the feed elements along the arms. 
By tracking the field over 12 h, a full synthesis image can 
be formed. The near-continuous UV coverage from 15 m to 
1.6 km results in good response to complex structure and low 
sidelobe levels. Technical specifications are shown in Tab le CD 
The MOST has b een described in detail bv lMillsl i|l98ll ) and 
iRobertsonl i|l99ll ). 

Since 1986 the MOST has observed a single field for a 
12 h synthesis almost every week night and often during the 
day on weekends. At the beginning and end of each 12 h 
synthesis, a set of up to 8 different calibrator sources is ob- 
served. Calibrator measurements are discarded if they do 
not pass a number of checks, and the remaining measure- 
ments are averaged to obtain gain and pointing solutions 
for the beginning and end of the observation. The gain and 
pointing solutions are linearly interpolated over the synthe- 
sis time between two calibr ator observations. The full list 
of calib rators is described by ICampbell- Wilson fc Hunsteadl 
i|l994bl ). Known variable calibrators were removed from the 
list following the analysis of IGaensler fc Hunsteadl ([2000). 

Our present analysis is performed on the final images 
proces sed according th e proce dure described bv lGreen et al.l 
l| 19991 ) and lBock et all (|l999l ). 



2.1 Observing modes 

The MOST is capable of observing in a number of different 
modes depending on the desired signal-to-noise ratio and 
field size, as summarised in Table [2] The MOST field of 
view is elliptical with the major axis aligned in the North- 
South direction and size a — aocsc \ 8\, where S is the decli- 
nation of the field centre and ao is the mode-specific minor 
axis size. The area of an image is approximately equal to 
7r(oo/2) 2 csc \Si\ and the spatial resolution of each image is 
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Table 1. Technical specifications of the MOST. The Declination 
range is given for a fully synthesised (12 h) image. The sensitivity 
and field sizes are mode-dependent and shown in Table [2] 



Parameter 



Value 



Centre Frequency 

Bandwidth 

Polarisation 

Declination range for full synthesis 
Synthesised beam 
Restoring beam 
Area per pixel 
Dynamic Range (typical) 



843 MHz 
3 MHz 

Right Hand Circular (IEEE) 
-30° to -90° 
43 X 43 esc \ 5\ arcsec 2 
45 X 45 esc 1 5 1 arcsec 2 
11 X 1 1 esc | (S | arcsec 2 
100:1 



the same for each mode at 45 X 45 esc \S\ arcsec. Examples of 
images observed in two different modes are shown in Fig. Q] 

All MOST images are primary-beam corrected, which 
results in increased noise towards the edges of images. The 
increased noise is most pronounced in the wide-field T im- 
ages, which are larger than the half-power point of the pri- 
mary beam. Therefore, we define an ellipse-shaped available 
area, parameterised by the minor axis shown in Table [21 
outside which we ignore any detections and measurements. 

We define an image as 'usable' if it passes the visual 
inspection and post-facto calibration procedure described in 
Section 13.41 
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Figure 3. Area of epoch vs. number of epochs for blin d ra- 
dio transient surveys. At the ext remes. iBower et al.l ll2007fl cov- 
ered 0.03 deg 2 with 844 epochs. [ Levinson et al .1 J2002l 7" covered 
5990 deg 2 with 2 ep ochs. ICarilli et al. [120031 ) covered a 0.002 deg 2 
with 5 epochs and ICroft et al] l|201ol l covered 690 deg 2 with 12 
epochs. This work, shown in red squares, spans some of the in- 
tervening parameter space. VAST-Wide and VAST-GP are the 
proposed wide-field and Galactic plane surveys of the Variables 
And Slow Transients (VAST) proj ect, planned for A SKAP, and 
the Pi GHz Sky Survey (PiGSS) l lBower et al.ll2010l) is planned 
for the Allen Telescope Array. 



2.2 Archive coverage 

The primary survey conducted with MOST was the Syd- 
ney U niversity Molonglo Sky Survey (SUMSS; lMauch et al] 
l2003h , which covers the sky south of —30° excluding the 
Galactic plane (|6| < 10°), and was carried out between mid- 
1997 and 2006. The MOST was also used to conduct two 
Molo nglo Galactic Plane Surveys: MGPS-1 (|Green et all 
1999), covering the Galactic p lane from 245° ^ I 355° 
and |6| < 1.5°, and MGPS-2 ([Murphy et al] 120071 ) cover- 
ing 245° < I < 365° and |6| < 10°. There has also been 
a wide range of directed observations and monitoring pro- 
grams. Here we analyse 7227 images, observed between 1986 
December 18 and 2008 August 28. 

A plot of the total and usable sky coverage of the archive 
is shown in Fig. [2] All of the sky south of —30° has been 
covered at least once due to the SUMSS survey, and the 
Galactic plane has been observed twice with the MGPS-1 
and MGPS-2 surveys. Some regions have been covered multi- 
ple times due to monitoring programs (e.g. the microquasar 
GX339— 4 has been observed 14 times) or deep imaging (e.g. 
the Hubble Deep Field South, which has 52 usable epochs). 
A number of images that were excluded from further analy- 
sis due to poor quality discovered during visual inspection, 
or difficulty in obtaining a post-facto calibration, resulting 
in the usable coverage being significantly less than the total 
coverage. 

For the purposes of comparing survey areas, we define 
the two-epoch equivalent coverage area as: 

N pos 

— epoch = ^ Ai ( -A/epoch, i 1), (1) 
i=l 

where N pos is the total number of independent directions on 



Table 3. Two-epoch coverage area (A two _ epoc i 1 ) for a range of 
recent and proposed blind radio variability surveys. VAST-Wide 
and VAST-GP are the wide-field and Galactic plane surveys of 
the Variables And Slow Transients (VAST) project, planned for 
the Australian Square Kilometre Array Pathfinder (ASKAP). 



Survey ^two— epoch 

(deg 2 ) 



Carilli et al. (20031 0.0008 

Bower et al. (2007) 28 

This Work 2776 

Levinson et al. (2002) 5990 

Croft et al. (2010) 6900 

PiGSS (proposed) (Bower et al. 2010) 10000 

VAST-GP (proposed) 38250 

VAST-Wide (proposed) 7290000 



the sky measured in the survey, Ai is the area of a particular 
direction of interest, and A^poch.i is the number of epochs a 
particular direction has been measured. 

A comparison of the two-epoch equivalent coverage area 
of the MOST archive with other blind radio variability sur- 
veys is shown in Fig. [3] and Table [3] The MOST archive cov- 
ers ~ 10 4 deg 2 once, and ~ 2 x 10 3 deg 2 twice. The MOST 
archive has a maximum of 52 epochs on a small (~ 1 deg 2 ) 
patch of sky. At the large-epoch extreme. lBower et al] (|2007T 1 
covered a very small (0.024 deg 2 ) patch of sky 844 times. 
Plots of the two-epoch equivalent coverage area versus decli- 
nation and Galactic latitude of the MOST archive are shown 
in Fig. H 

The interval between epochs covers a wide range. Fig- 
ure [5] shows a histogram of the number of intervals between 
non-redundant pairs of measurements of each source as a 
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Figure 1. Examples of MOST images from the MGPS surveys. Left: An T field from MGPS-2, Right: A 'B' field from MGPS-1. Both 
images are centred on I = 333.9° , b = —4.4° . The grey scale is the same for each image ranging from —5 mjy beam" 1 to +15 mjy beam - 1 . 
The resolution of each image is the same (45 X 61) arcsec, and the measured background RMS values are 1.5 and 0.85 mjy beam -1 
respectively. The 'B' field overlaps the T field in the region shown. A number of typical artefacts are visible, including radial spokes from 
bright sources, regions of negative emission, grating rings and increased noise towards the edge of the T field. 



Table 2. MOST observing mode, image parameters (see Section [27TJ and number of images in the archive. The columns are: the mode 
code, typical Root Mean Square (RMS) of the noise in the image, full and available minor axes, total number of images in the archive and 
equivalent area (assuming the available minor axis) and number and equivalent area of images for which visual inspection and post-facto 
calibration were s uccessful (see Section |3.4| |. The 'H' and T modes were only available since the 1993 wide-field upgrade described by 
lLarge et al] Jl994h . 
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Count 


Area 


Count 
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0.0 
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2.0 


1.36 


1.02 


3915 


18156.4 
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9266.1 



function of the interval in days. At the two extremes, Fig- 
ure [S] shows that there are over two thousand pairs of mea- 
surements separated by one day, and the same number of 
measurement pairs separated by over 20 years. The archive 
is most sensitive to changes in flux density on the time-scale 
of around one thousand days, where it contains over ten 
thousand pairs of measurements. The results for transient 
and variable radio sources are discussed in Section T5, 21 



2.3 Image artefacts and systematic errors 

2.3.1 Visually obvious artefacts 

MOST images suffer from a number of image artefacts that 
make accurat e automatic point source detection and fit- 
ting difficult (|Mauch et al.ll2003l ). In particular, some im- 
ages are contaminated by extended emission (especially in 
the Galactic plane), with negative bowls or radial spokes 
around bright sources, grating rings due to the telescope ge- 
ometry and bands across images due to solar and terrestrial 
interference (see, for example Fig. [TJ. To mitigate the ef- 
fects of these artefacts, we discard clearly affected images 
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Figure 2. Coverage of the MOST archive in a Lambert equal area projection sampled on 0.1°grid and centred on the South Celestial 
Pole. Left: the entire archive. Right: usable images after visual inspection and post-facto calibration (see Section l3.4l l. Some 4200 of the 
7200 images were not useable and appear only in the left panel. The colour axis is number of epochs with light grey being no images, 
dark grey being one image (unusable for detecting transients and variability) and colours yellow to blue indicating 2 or more images. 
The shape of each image is the available size shown in TablefS] not the full image size. The SUMSS survey covers the sky south of —30° 
and the Galactic plane surveys: MGPS-1 and MGPS-2 can be clearly seen as a band from the left-middle to top-right of both images. 
Repeated fields include the Hubble Deep Field South (52 usable epochs), SN1987A (13 usable epochs), GX339-4 (13 usable epochs) and 
GRO1655-40 (15 usable epochs). 



by visual inspection, considered only point source fits with 
small formal errors, and measured the RMS and mean in 
the region around each detection to correct for localised off- 
sets and noise. Extended emission was not explicitly nagged 
but point sources with small formal errors of fit were rarely 
found close to regions of extended emission. 

2.3.2 CLEAN artefacts 

As with any interferometer, poor dirty beam modelling can 
result in negative flux around sources after CLEANing (see 
I Green et al .11 19991 for a discussion). An example of this ef- 
fect is shown in Fig. [6] When performing automatic source 
fitting, the negative surface brightness around the edge of 
the source reduces the integrated flux density to a value 
significantly less than the peak flux density, which should 
not occur if CLEAN is using the correct beam. In a sam- 
ple of 6767 otherwise good point source fits in the range 
h < RA < 10 h, only 45 per cent had peak and integrated 
flux density measurements in formal agreement, and more 
than half of these had a peak value that significantly ex- 
ceeded the integrated flux density. Assuming that variable 
and transient sources must be point sources at the MOST 
resolution, integrated flux densities were not used in the suc- 
ceeding analysis, with peak values used throughout. 

2.3.3 Variable calibrators 

iGaensler fc Hunsteadl ([2000) found that, between 1984 and 
1996, 18 of the 55 Molonglo calibrators were variable on 
characteristic time-scales between 300 and 3000 days. The 
Molonglo calibration procedure averages the measurements 
of up to 8 calibrator sources, so variation of a single source 



is expected to have a relatively small effect on the flux den- 
sity scale calibration. The characteristic time-scale was much 
longer than the synthesis time so we expect no variability on 
that time-scale. In 1996, the variable sources were removed 
from the calibrator list and the post- 1996 flux density cali- 
bration (which included SUMSS and MGPS-2) can be con- 
sidered more reliable. In any case, the flux density scale for 
every image was calibrated against SUMSS/MGPS-2 using 
the post-facto calibration procedure (see £|3.4fl , so any resid- 
ual error in the flux density scale should be removed by this 
technique. 

2.3.4 Radio frequency interference 

Observations performed after 2007 January at the Molon- 
glo site have been subject to interference from strong local 
sources of Radio Frequency Interference (RFI). The inter- 
ference affects almost all wide- field T images observed after 
2007 January, but 'B' fields are much less affected. No other 
modes were used during this period. The strongest RFI is 
usually excised during the reduction process, which results 
in incomplete coverage of the UV plane and higher sidelobe 
levels. Remaining low-level RFI resulted in significantly in- 
creased image noise levels. Occasionally, the interference can 
manifest itself as gain variations across the image, which is 
difficult to identify with visual inspection alone. Images af- 
fected in this way were detected and removed by the post- 
facto calibration procedure (see ij3.4|) . 

2.3.5 Missing calibration values 

The standard MOST image reduction pipeline assumes de- 
fault values for calibration parameters if none are supplied. 
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Figure 4. Two epoch coverage vs. declination (above) and vs. 
Galactic latitude (below). The horizontal error bars denote the 
bin widths of 10°. 
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Figure 5. Histogram of the number of non-redundant pairs of 
measurements in the light curves of all sources in the sample, the 
variable sources only, and the transient sources only, binned by 
the interval between the measurements in days. 
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Figure 6. Example of poor dirty beam modelling producing neg- 
ative regions around SUMSS J022129-335825 after CLEANing. 
Fitting a 2D Gaussian in a 10x10 pixel region centred on the 
source results in an integrated flux density of 394.0 mjy which is 
less than peak intensity of 449.3 ± 16.6 mjy beam -1 . The size of 
the fitted Gaussian is consistent with the size of the synthcsiscd 
beam (shown in the bottom left). The large difference in inte- 
grated and peak flux densities is due to poor dirty beam mod- 
elling. 



Inspection of image headers revealed that over 600 images 
had been reduced with default calibration values, rather 
than the measured values. We were able to re-reduce these 
images with the sa me calibration measurements used to pro- 
duce the results of iGaensler fc Hunsteadl (|2000l ). 

In spite of this effort, we found large variations in the 
flux density scale on some occasions, which we attribute to 
the variable calibrators, wrongly applied calibration and in- 
terference during calibrator observations. Rather than dis- 
carding images that appeared good apart from flux density 
scale, we included them in our analysis and calibrated with 
the post-facto calibration procedure (see i]3.4fl . 

2.3.6 Resolution effects 

Many radio t ransient and va r iabili ty surveys, such as that 
performed bv lLevinson et al.l(|2002l ). have suffered from the 
potential for mismatched resolutions to cause systematic er- 
rors in the analysis. Mismatched resolutions can cause a 
number of systematic effects. In particular, higher resolu- 
tions can resolve a large source into multiple components, 
as well as resolve out the extended flux, resulting in a lower 
flux density measurement. 

The MOST synthesised beam is fixed in the east-west 
direction, but elongates in the north-south direction with the 
declination of the image centre. The difference in resolution 
is largest when a source is imaged at the southern tip of an 
image centred on the declination limit (f5 max = —30°), and 
the northern edge of an image centred to the south, with dec- 
lination S2 — <5max — 2ao,max esc |5 ma x|- The largest minor axis 
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occurs for the T images, with. <2o,max — 1.36° (see Table©. 
The difference in resolution in the worst case is therefore 
45(csc |5max| — esc |($2|) ~ 9.2 arcsec in the north-south di- 
rection, which corresponds to a maximum difference 
in resolution of overlapping images of approximately 

20 per cent. 

By way of comparison. iLevinson et al.l lj2002T l compared 
NVSS and FIRST surveys, whose resolutions differed by a 
factor of 5 (45 to 5 arcsec). We suspect resolution difference 
affects a very small fraction of our measurements as only 
those sources in the narrow band of overlap between north- 
south adjacent images at northern declinations are affected, 
and that multiple match and extended flux effects will be 
small as the change in resolution is only 20 per cent. 



3 LIGHT CURVE EXTRACTION 

There are two broad categories of techniques that can be 
applied to a set of astronomical images to detect vari- 
able and transient emission: image comparison and cata- 
log ue comparison. Im age comparison, as used for example 
bv lCarilli etafl (|2003h . involves comparing an image of in- 
terest with a reference image, either by subtraction or di- 
vision. Source finding can then be performed on the differ- 
ence or quotient image. Catalogue comparison makes source 
catalogues by performing source finding on each image and 
forming a light curve by measuring the flux density sepa- 
rately on each imag e. The light curve i s then subjected to 
a statistical analysis. iBower et al l 1120071 ) used a cataloguing 
approach. 

An image comparison approach is not feasible for the 
MOST archive since MOST image artefacts change as a 
function of pointing centre and observing mode, which would 
result in a large number of false detections when a reference 
image is compared with a candidate image. Instead we use a 
catalogue-based approach as discussed in the following sec- 
tions. 



3.2 Source association 

An initial catalogue was created by merging the SUMSS 
catalogu^U and the MGPS-2 catalogu^El an d source detection 
was performed independently on each image as described 
in § 13.11 We associated each detection with a catalogued 
source if the positions of the detected and catalogued sources 
were coincident within 10". If there was no previously known 
source within 10", a new source was added to the catalogue. 

3.3 Measurement 

Once a source has been detected, we measured flux density 
at the position of the source on all images in the archive for 
which this position is in the field of view. We do this for a 
number of reasons: 

(i) To obtain fitted flux densities at all epochs. 

(ii) To a obtain upper limits on flux density for epochs in 
which a source is not present. 

(iii) To obtain measurements of the surface brightness in 
a region around a source, which is useful for quantifying 
the RMS error of the measurement (for error estimates) and 
correcting for a non-zero mean level. 

To obtain a post-facto calibration solution for an image, 
we require accurate flux densities of bright point sources. 
To measure flux density for bright sources we fit a Gaus- 
sian source to a 10 x 10 pixel square box on the detected 
position. For weak sources (< 5a), Gaussian source fitting 
occasionally produces large errors, partly due to the large 
number of free parameters. Therefore, to measure flux den- 
sity of weaker sources, we also performed a 2D parabolic 
fit to a 3 x 3 pixel region centred on the source position. 
To quantify a non-detection, we measured the peak pixel 
in a 3 x 3 pixel square region centred on the source posi- 
tion, and to measure the local background we measured the 
RMS and mean flux density in an annular region of inside 
radius 66" and outside radius 99" (6 and 9 pixels in radius, 
respectively). We did not mask out other bright sources in 
the annular region, so measurements of the local background 
can be biased in densely populated areas. 



3.4 Post-facto calibration 



3.1 Source detection 

Source detection w as performed using the SFIND 

task (jHopkins et al.l I2002T ) from the miriad package 
|Sault fc Killeenl I2009T ). SFIND uses the False Discovery 
Rate method for controlling the fraction of false positives 
when performing multiple hypothesis testing. It accepts two 
parameters: the size of the square box over which it should 
measure the background RMS surface brightness, and the 
percentage of false pixels to accept when applying the False 
Discovery Method. We chose a box size of 50 x 50 pixels 
and a false pixel rate of 10 per cent. SFIND automatically 
produces Gaussian source fits and background RMS for 
each detected source. We used only the fitted positions. We 
separately measured the flux densities as described in N3.3I 
below as we found the fitted flux densities and background 
RMS from SFIND to be occasionally unreliable. 



In our initial investigations, we found a large number of 
sources that appeared to have varying flux densities, which 
we found was due to variations in the gain of individual 
images. 

In order to reduce the number of false variability de- 
tections due to this effect, we implemented a post-facto 
calibration scheme by comparing the Gaussian-fitted peak 
flux density measurements of the bright sources on each 
image, against th e SUMSS flux densities measured by 
iMauch et ail (120031) or t he MGPS-2 flux densities measured 
bv lMurphv et alj (|2007f l. 

For each image, we performed a 2-stage fit for the gain. 



1 Version 2.1 dated 2008 Mar 11 

http : //www. physics .usyd. edu. au/sif a/Main/SUMSS 

z Version dated 2007 Aug 15 

http : //www. physics .usyd. edu. au/sif a/Main/MGPS2 
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Table 4. Threshold values of fitted formal errors used to select 
point sources for detection and post-facto calibration. 



Parameter 



Formal error threshold 



Major axis 
Minor axis 
Position in RA 
Position in Dec 
Peak Value 



< 15" 

< 12" 
<5" 
< 6" 

< lOmJy beam -1 



First, we selected all point sources that had formal fit errors 
within the bounds shown in Table [4] a peak flux density 
above 40 mjy, and a counterpart in SUMSS or MGPS-2. We 
then performed a linear least squares fit to the relationship 
between peak fitted flux density measurements and cata- 
logued flux densities, weighted by the inverse of the quoted 
peak flux density measurement error (a) from the catalogue. 
The fit was constrained to go through zero. We use the fit- 
ted Gaussian measurements because it provides a more ac- 
curate flux density measurement than the 2D parabolic fit 
for the bright sources that we used for the calibration and is 
more sensitive to image errors, thus enabling the rejection 
of images containing subtle image errors that affect a large 
fraction of sources. All points that deviated from the fit by 
more than 3a were removed from the set and the fit was per- 
formed again on the remaining points. The gradient of this 
fit was used as the gain calibration factor, and was applied 
to all measurements derived from the image in question. 

In addition to producing a gain estimate, this proce- 
dure also enables us to remove poor quality images from the 
remaining analysis by measuring the goodness-of-fit, both 
before and after the deviant points were removed. As a mea- 
sure of goodness-of-fit we adopted the reduced % 2 1 defined 
as 



-2 _ Xce.\ 
Acal 



Xcal 



N d 



iVcal - 1 



(2) 



where Adof is the number of degrees of freedom in the fit and 
Acai is the number of sources being used for the calibration. 

We discarded images for which the following conditions 
were all met: 

• The initial Xcal was greater than 30. 

• More than 30 per cent of the sources were more than 
3cr from the fit. 

• The x 2 was greater than 5 after removing the 
sources which were 3a from the fit. 

The post-facto calibration procedure assumes that the 
majority of sources in the sky are static over time-scales of 
several years and variations in flux density are not spatially 
correlated over the MOST field of view. As such, refractive 
scintillation caused by a nearby cloud of ionised gas cover- 
ing a significant fraction of the MOST field of view will be 
removed by the calibration and rendered undetectable. We 
consider this possibility physically unlikely. 

Examples of the post-facto calibration procedure ap- 
plied to two images are shown in Fig. [7| A plot of the Xcal 
vs time is shown in Fig. [9] from which we were able to de- 
tect the interference that began affecting the telescope in 
2007 by a group of images with large Xcai- This increased 
our confidence in the usefulness of the post-facto calibration 
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Figure 7. Examples of post-facto calibration of two images. Cat- 
alogued flux densities (from SUMSS/MGPS-2) are plotted against 
the peak flux density measured with a Gaussian fit on a 10x10 
pixel box centred on each source position. The solid line is the 
line of best fit after removing the sources more than 3<r from the 
initial best fit line including all points (where a is the quoted er- 
ror from SUMSS/MGPS-2). The squares arc the points included 
in the fit and the crosses are the sources removed. Top panel: 
An T image (observed 2000 March 25) with good fit and gain 
of 0.89. Bottom panel: An T image (observed 1999 November 9) 
with a poor fit due to large discrepancies between measured and 
catalogued flux densities for many sources, and gain of 0.98. 



procedure and, in particular, in its ability to recognise poor 
quality images. 

An example of the post-facto calibration procedure ap- 
plied to the light curve of a single source is shown in Fig. [8] 
The procedure removes the most discrepant points and cor- 
rects the remainder, which results in less scatter in the light 
curves. 



3.5 Light curve characterisation 

3.5.1 Variability Classification 

Astronomical sources exhibit a wide variety of light curve 
shapes. For example, interstellar scintillation produces ran- 
dom variations around a mean, whereas stellar flares are 
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Figure 8. Light curves illustrating the effect of post-facto cal- 
ibration on SUMSS J201139-561858. Top panel: Pre-calibrated 
light curve, centre panel: calibrated light curve, bottom panel: 
difference between the calibrated and uncalibrated light curves. 
The uncalibrated light curve has a modulation index (see Equa- 
tion |9j of 0.181 compared to that of the calibrated light curve of 
0.0842. Two images were discarded from the light curve by the 
post-facto calibration. 



episodic with short duty cycles and flux densities signifi- 
cantly above a quiescent level. 

Determining whether a given light curve matches any 
of a wide variety of light curve shapes, and quantifying 
the fidelity of the match, are beyond the scope of this 
work. For simp l icity, we adopt the classification measure of 
iKesteven et al.l (1 19771 ) to measure random variability in the 
light curve, namely the x °f the residuals of a weighted fit 
to a line of constant flux density. It is computed according 
to the following expressions: 



Xlc 



E 



(Si - By 



where: 
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Figure 9. x 2 , of the post-facto calibration fit versus image ob- 
servation date, before and after 3<r outliers were removed. Note 
the significant number of images with large value of reduced Xcni 
(before outlier removal) in 2007 and 2008, when terrestrial inter- 
ference was affecting the telescope. 
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Figure 10. Histogram of post-facto calibrated gains for 2 999 
usable images. The dashed line is a Gaussian distribution fit to 
the gains, with a mean gain of 0.985 and standard deviation of 
0.056. 



/ rirav 

= g{Si 



Pi) 



(6) 



and S is the weighted mean flux density, <7; is the 
estimated flux density error, Sl aw is the flux density 
measured on the image with the 2D parabolic fit, g is the 
gain coefficient calculated by the post-facto gain calibration 
and pi and <7b g ,i are respectively the mean and RMS sur- 
face brightness of an annular region around the source. The 
weighting factor of 1 /erf captures both the noise level in the 
surrounding region and the absolute flux density calibration 
error, from the MGPS-2 and SUMSS catalogues, assumed to 
be 5 per cent. For a source with a peak flux density at the de- 
tection limit of 14 mjy beam -1 , ai is dominated by the back- 
ground RMS term and has a typical value of ai ~ 3mJy. 
For a source with flux density of 100 mjy beam -1 the flux 
density calibration term dominates and results in a typical 
error of ai ~ 6mJy. 



10 K. Bannister et al. 



To determine whether a source is variable or not, we 
work from the null hypothesis that all sources are static 
with added Gaussian measurement error. Applying the mea- 
sured Xic t° t ne analytic x 2 cumulative distribution function 
(CDF) yields the probability of rejecting the null hypothe- 
sis, i.e. the probability that the source is genuinely variable, 
which we write as P(xt c )- Working with probabilities rather 
than x 2 values also enables us to compare the statistics of 
light curves with differing number of measurements, which is 
absorbed in the CDF. Those light curves whose probability 
exceeds a certain threshold can be classified as variable. 

Pa st investigations such a s th ose of Kesteven et al.l 

l| 19771 ) .IC rcgorv fc Tavlorl ij 19861 ) and lSeielstad et al.l (|l983T ) 
applied a probability threshold of 99 per cent to derive a se t 
of variable sources. In the case o f iGregorv fc Tavlorl ljl986h . 
the 99 per cent threshold left 4 per cent of the sources clas- 
sified as variable. Of these sources, 1 per cent of the total, or 
25 per cent of the classified variables must be false positives, 
simply from the choice of a low threshold. 

As we have almost 30 000 light curves, a 25 per cent 
false positive rate is too high. To arrive at a more rigourous 
threshold we define a confidence measure that is the excess 
sources over expected number given the null hypothesis, ex- 
pressed as a percentage. Specifically, we use: 



C{T) 



N n . 



i(T) — N a 



aly(T) 



(7) 



where Actual (T) is 

,2 \ 



the number of sources with light 
curves whose -P(x 2 c ) is greater than T, Anoisconiy (T) = 
A^ourccs (l — P (xfc > ^)) i s the number of sources expected 
if there were only noise in the sample, and Asourccs is the to- 
tal number of sources in the sample. The confidence measure 
C(T) in Equation [7] represents the fraction of genuinely vari- 
able sources contained in the set of sources with P(X\ C ) > T. 



3.5.3 Detections 

A measurement must satisfy a number of criteria in order to 
be considered an acceptable measurement of a point source 

rather than a result of an image artefact. 

The analysis of the SUMSS catalogue bv lMauch et al.l 

(2003) showed that the SUMSS survey was complete above 
15 mjy and we found a number of convincing detections 
above 14 mjy. Gaussian fits can have major and minor axes 
significantly different from the synthesised beam, either be- 
cause a source is partially resolved, or because of a strong 
image artefact. Restricting the major and minor axes of the 
fit to less than 20 per cent larger than the synthesised beam 
ensures the fit is a point source and reduces the number of 
false detections due to artefacts. 

Therefore, we define a 'detection' as a Gaussian fit to a 
source that satisfies all of the following criteria: 

• Fit errors are below the thresholds in Table [3] 

• Fit position is less than 5" from the catalogued position 

(c.f. §EU), 

• Peak flux density ^ 14 mjy, 

• Ratio of fit major axis to synthesised beam major axis 
is < 1.2, 

• Ratio of fit minor axis to synthesised beam minor axis 
is < 1.2, 

• Image passes the post-facto calibration. 



3.5.4 Detection signal-to-noise ratio 

To find transient sources, it is particularly important to 
measure the significance of a given detection or non- 
detection. We define the detection signal-to-noise ratio as 
d = Sccnmax/o-tg where Sccnmax is the maximum surface 
brightness in a 3 x 3 pixel box centred on the source po- 
sition, and (Tb g is the RMS surface brightness (in mjy) in 
an elliptical region centred on the source position. 



3.5.2 Variability measures 

Previous work has adopted a number of different metrics 
to quantify variability. To enable comparison with previous 
work, we adopt two different metrics, the first is the frac- 
tional variability defined as: 

V*Smax 5mm . * 

- -5 T O . ■ W 

'-'max r '-'mm 

where S max and Smin are respectively the maximum and 
minimum flux densities. This measure is most suited to 
characterising sources wi th few measurements or f laring be- 
haviour and was used bv lGreeorv fc Tavlorl l|l986h . 

The second metric is the modulation index, defined as: 

m = a/B (9) 

where a is the standard deviation of the light curve and 
S is the weighted mean defined in Equation [4] This met- 
ric is more suited to measuring random variability in data 
sets with a large num b er of measurements and was used by 
iGaensler fc Hunsteadl l|2000l ). 



3.6 List of variable and transient sources 

We computed all metrics described in M3.5I for all sources 
and produced lists of variable sources, and another list of 
transient sources. The list of variable sources was produced 
by finding all sources with: 

• Light curve probability exceeding a threshold that 
yields a confidence metric greater than 90 per cent, i.e. 
P{xl) > To where C{T ) > 90 percent. 

• A 'detection' on every image. 

A variable source, therefore, has at least two detections 
and no non-detections. The list of transient sources was pro- 
duced by finding all sources with: 

• Minimum detection signal-to-noise ratio d m in < 5.0 

• Maximum detection signal-to-noise ratio d max > 6.0 

• Flux density at the maximum detection signal-to- 
noise ratio epoch (measured with a parabolic fit) > 
14 mjy beam -1 

• P(xic) > 99% computed with Scenmax as the flux den- 
sity value. 

A transient source, therefore, has at least one detection 
and at least one non-detection. We use the parabolic fit, 
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Figure 11. Histogram of the number of epochs measured per 
source in the variability and transient analysis. At the two ex- 
tremes, there are 26000 sources with two measurements 
and 11 sources with 51 measurements. Gaps in the his- 
togram indicate no sources have that number of epochs. 



as this was less prone to error at low signal-to-noise ratio. 
We compute the Xic f° r the transient list using the Scenmax 
rather than the 2D parabolic fitted flux density as this is the 
most robust measure of the upper limit to the flux density, 
and we apply the error as per Equation [5] We reduce the 
threshold on the P(xfc) with respect to the variable sources, 
because the flux density of the non-detections is not known. 



4 RESULTS 

Throughout this section we assume a Hubble constant of 
Ho = 72 km s _1 Mpc _1 and define the spectral index, a, 
according to S oc v a where S is the flux density and v is the 
observed frequency. 

When referring to a radio source, we use the popular 
name (e.g. SN1987A) if one exists. Otherwise, if the source 
appears in either the MGPS2 or SUMSS catalogues we use 
the MGPS2 or SUMSS naming convention, or if not, we use 
the source coordinates in the form 'JHHMMSS-DDMMSS' 
without prefix. 



4.1 Data Quality 

4-1.1 Number of usable images 

Of the 7227 images in the archive, 2999 were retained after 
visual inspection and post-facto calibration. The area of an 
equivalent two-epoch survey (see Equation[T]) is 2775.7 deg 2 . 



4-1.2 Number of measurements 

There are 29730 sources that have light curves with measure- 
ments at two or more epochs. A histogram of the number of 
epochs measured on each source is shown in Fig. 1111 



4-1.3 Analysis of the population of light curves 

To check that the flux density error estimate is consistent 
over the population of light curves, we compare the mea- 
sured and theoretical distributions of the residuals from a 
line of constant flux density. As described in Equation [3] 
the residuals of the best-fitting light curve of constant flux 
density follow a x 2 distribution with iV m eas — 1 degrees of 
freedom. Fig. ll2l shows the histograms of the x? c of the resid- 
uals for light curves in our variability sample. The theoreti- 
cal and actual curves agree reasonably well, with the slightly 
larger number of sources at small values of xfc indicating 
that the flux density errors have likely been marginally over- 
estimated. 

In order to compare distributions of sources with differ- 
ing numbers of measurements, we can compute the proba- 
bility of a given xfc value occurring by chance by using the 
analytic \ 2 cumulative distribution function with the appro- 
priate degrees of freedom. If the variation in the light curves 
is attributed entirely to Gaussian noise, the histogram of the 
computed probabilities should be uniform and independent 
of the degrees of freedom. Fig. fTS] depicts the expected and 
actual probability histograms. There is a general trend in 
the actual histogram showing an excess of sources with high 
1 — P(Xlc)> which we attribute to overestimated flux den- 
sity errors leading to smaller measured xL values. Because 
a large fraction of sources have a small number of measure- 
ments (see Fig. Hip , at this degree of freedom a small Xic 
corresponds to a higher value of 1 — P(xic) which in turn 
leads to the observed excess. In addition to the general trend, 
there is a single bin at the lowest 1 — P (xi c )> which we at- 
tribute to a population of genuinely variable sources. The 
genuinely variable sources are clearly seen as an excess of 
sources, departing from the theoretical curve as the proba- 
bility is reduced. 

Fig. ll4l shows the confidence measure C(T) vs 1— P(xi c ) 
on a log scale. From this plot we deduce that applying a 
threshold of T = 1 - P(xic) = 4 X W 5 

ensures the fraction 

of spurious variable sources in the list of candidate variables 
is less than 10 per cent. 



4-1-4 Primary beam effects 

In the MOST archive, a significant fraction of the over- 
lapping coverage is at the edges of adjacent images. Er- 
rors in primary beam correction, which increase with ra- 
dial distance from the pointing centre, could therefore in- 
troduce false variability. To investigate this effect, we took 
the light curves of all sources appearing in images centred 
on a =22h30m00s, S =-61d06m36s (J2000), in the Hubble 
Deep Field South (HDFS) region. We calculated the fraction 
of measurements that deviated by more than 20 per cent 
from S (see Equation binned by radial distance from the 
pointing centre, as shown in Fig. 1151 We found no strong 
dependence with distance and therefore conclude that the 
primary beam correction is accurate, and that the adopted 
available minor axis of the primary beam (see Tableland § 
I2,l[l . outside which we ignore measurements, is acceptable. 
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Figure 12. Histograms of \f c values for light curves with between 2 and 10 measurements (Nmeas)- Also plotted as filled circles is the 
theoretical v 2 ,, , probability distribution with A^of = Wmcas — 1 degrees of freedom binned as for the measurements. 



Transient and variable radio sources with MOST 13 



1 




■ Actual 
Expected 






V 



l-f(Xl 2 c 




Figure 13. Top panel: Number of sources vs 1 — P(xf ) on a 
linear scale. The y error bars are the Poisson errors and the x error 
bars indicate the bin widths. The solid line depicts the expected 
number of sources in each probability bin if the scatter in the 
light curves was due to Gaussian noise. The point at the far left 
indicates a population of variable sources and the excess over the 
solid line indicates that there are approximately 350 genuinely 
variable sources in our sample. The increasing general trend from 
left to right indicates there is a deficiency of moderately variable 
and excess of static sources, which agrees with Fig. 1121 The region 
to the left of the dashed vertical line at 1 — P(xf c ) = 0.2 indicates 
the baseline number of sources used to infer the actual number of 
variable sources (see i|5.5ll . Bottom panel: The same data but 
re-binned on a log scale to illustrate how the number 
of sources departs from the expected number for small 
values of 1 — P(x? )• 
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Figure 14. Confidence measure (C(T)) vs 1 — P(x? ) on a log 
scale. To achieve a fraction of genuine variable sources of 90 per 
cent requires 1 — P(xi c ) = 4 X 10 — 5 . The y errors are the Pois- 
son errors propagated through the definition of the confidence 
metric, and the the x errors indicate the bin width. The confi- 
dence is negative at large values of 1 — P(x? c ) due to a deficiency 
of moderately variable sources at large values of 1 — P(Xi 2 c ) (see 
Fig. 1131 ). The dashed vertical line indicates the adopted threshold 
of 1 — P(xf c ) = 4 X 10~ 5 = To, to the left of which the confidence 
of finding a genuinely variable source is greater than 90 per cent. 
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Figure 15. Plot of the fraction of flux density measurements that 
deviate from S by more than 20 per cent, binned by distance 
from the image centre. Only images centred on the HDFS were 
considered, and 2 known variable sources were removed. The error 
bars are Poisson errors 



4- 1.5 Declination effects 

The sensitivity and resolution of MOST images is strongly 
declination dependent, with noise and beam size increas- 
ing towards Northerly declinations. Therefore, to ensure our 
sample is not contaminated by declination dependent sys- 
tematic effects, we plot in Fig. \TE\ the fraction of transient 
and variable sources as a function of declination. We find no 
declination dependence. 



4.2 Variable sources 

Table [5] lists the 53 sources that we classify as variable ac- 
cording to the criteria described in H3.6I The Table also lists 
radio and optical counterparts where available. Light curves 
for each of these sources are shown in Appendix 1X1 
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Figure 16. Fraction of variable and transient sources vs declina- 
tion. 



Table 5. 53 radio sources classified as variable with a confidence greater than 90%. The columns are: the source name, right ascension and declination (J2000), Galactic longitude and 
latitude, number of accepted measurements, probability of \\ °f the residuals to a line of constant flux density, maximum, weighted mean and minimum flux density at 843 MHz, 
modulation index, V value, NVSS flux density at 1.4 GHz (if covered by NVSS area), SuperCOSMOS classification: D = point source detected, N = no source detected, C = more than 
1 point source in the error box, E = extended according to SuperCOSMOS catalogue, G = clearly resolved galaxy, the SuperCOSMOS optical magnitudes in blue (B), red (R) and 
infra-red (I), and the B— R colour. 
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20.46 


— 




.48 


SUMSS 


J003722— 421508 


00 


37 


22 


-42 


15 


12 


312 


8 


-74. 


6 


2 


6.5 X 10~ 6 


25.3 


19 


.2 


15.2 





.4 





.25 


- 


D 


20.54 


20.15 


— 




.38 


SUMSS 


J004417-375258 


00 


44 


17 


-37 


52 


57 


310 


4 


-79 


1 


2 


2.4 X 10 -14 


43.5 


23 


.8 


18.7 





.6 





40 


18.1 


E 


20.30 


19.14 


18.91 


1 


.16 


SUMSS 


J010426-400412 


01 


04 


26 


-40 


04 


12 


292 





-76 


8 


2 


2.3 X 10 _s 


62.6 


44 


.6 


37.7 





.4 





25 


- 


D 


20.00 


19.15 


18.92 




.85 


SUMSS 


J011019— 455112 


01 


10 


20 


-45 


51 


16 


292 


8 


-70. 


9 


7 


2.3 X 10~ 6 


20.9 


16 


.0 


10.7 





.2 





.32 


- 


G 


— 


— 


— 




— 


SUMSS 


J011136— 394340 


01 


11 


36 


-39 


43 


40 


285 


8 


-76 


7 


2 


1.0 x 10 -9 


44.5 


27 


.9 


23.5 





.4 





.31 


21.5 


E 


19.32 


18.35 


18.77 




.97 


SUMSS 


J011839— 452657 


01 


18 


40 


-45 


26 


59 


288 


2 


-70 


9 


9 


9.4 x 10~ 6 


45.7 


37 


4 


32.0 





1 





.18 


- 


N 


— 


— 


— 




— 


SUMSS 


J012615— 561946 


01 


26 


15 


-56 


19 


47 


293 


2 


-60 


1 


2 


9.5 X 10 -6 


14.1 


9 


.3 


6.4 





.5 





.38 


- 


C 


— 


20.50 


18.35 




— 


SUMSS 


J013013 — 742024 


01 


30 


15 


-74 


20 


24 


299 


4 


-42. 


5 


4 


1.0 x 10 -5 


29.9 


24 


.2 


20.9 





2 





.18 


- 


N 


— 


— 


— 




— 


SUMSS 


J021036— 405154 


02 


10 


36 


-40 


51 


55 


258 


1 


-68 


7 


2 


8.8 X 10~ 6 


40.4 


32 


.6 


26.1 





.3 





.22 


- 


D 


21.28 


19.96 


— 


1 


.32 


SUMSS 


J022223— 802356 


02 


22 


23 


-80 


23 


59 


298 


4 


-35 


9 


2 


2.4 x 10~ 7 


56.7 


41 


1 


35.0 





.3 





.24 


- 


D 


19.95 


19.70 


— 




24 


SUMSS 


J031201-523432 


03 


12 


00 


-52 


34 


32 


266 


9 


-53 


6 


2 


1.4 x 10~ 5 


30.8 


24 


.7 


18.9 





.3 





.24 


- 


N 


— 


— 


— 




— 


SUMSS 


J032102— 405735 


03 


21 


02 


-40 


57 


32 


247. 





-56 


4 


2 


2.9 x 10~ 8 


39.8 


27 


.8 


19.4 





.5 





.35 


- 


D 


20.21 


19.64 


— 




.58 


SUMSS 


J040608- 540445 


04 


06 


08 


-54 


04 


54 


263 


9 


-45. 


6 


2 


1.6 x 10~ 8 


19.6 


12 


.5 


9.0 





.5 





.37 


- 


C 


21.71 


19.33 


18.49 


2 


.38 


SUMSS 


J043624-563332 


04 


36 


24 


-56 


33 


31 


265 


7 


-40 


8 


3 


9.8 x 10~ 10 


81.4 


61 


.5 


49.6 





.2 





24 


- 


D 


20.22 


19.63 


— 




.58 


SUMSS 


J044508-541155 


04 


45 


08 


-54 


11 


53 


262 


4 


-40 





2 


1.3 X 10 _s 


45.7 


33 


1 


28.2 





.3 





.24 


- 


D 


20.71 


19.99 


— 




.72 


SUMSS 


J044624— 391548 


04 


46 


24 


-39 


15 


45 


242 


8 


-40 


3 


2 


2.5 X 10 -5 


47.1 


37 


.8 


32.2 





3 





19 


44.2 


D 


18.41 


17.64 


17.21 




.77 


SUMSS 


J044829-401714 


04 


48 


28 


-40 


17 


12 


244 


2 


-39 


9 


2 


3.4 x 10~ 9 


58.9 


43 


.8 


35.4 





.4 





25 


- 


D 


19.44 


19.39 


19.15 




.05 


SUMSS 


J045530— 720442 


04 


55 


30 


-72 


04 


41 


283 


9 


-34. 


5 


3 


< 1 x 10~ 14 


38.5 


17 


.0 


11.9 





6 





53 


- 


C 


21.50 


21.68 


— 


- 


.18 


SUMSS 


J050230-651321 


05 


02 


30 


-65 


13 


22 


275 


6 


-35 


8 


4 


2.3 X 10 -9 


40.4 


31 


.4 


25.5 





2 





.23 


- 


c 


19.47 


18.42 


— 


1 


.05 


SUMSS 


J061952— 643917 


06 


19 


52 


-64 


39 


17 


274. 


3 


-27 


7 


4 


2.3 X 10 _s 


49.5 


33 


.5 


28.5 





.3 





27 


- 


E 


21.57 


19.41 


19.32 


2 


.16 


SUMSS 


J062724— 405328 


06 


27 


24 


-40 


53 


28 


248 


9 


-21 


6 


2 


2.6 X 10~ 5 


24.0 


16 


.9 


13.7 





.4 





27 


- 


N 


— 


— 


— 




— 


SUMSS 


J064951— 620221 


06 


49 


51 


-62 


02 


21 


272 


1 


-23 


9 


2 


2.4 x 10~ 6 


46.2 


36 


.4 


30.1 





.3 





21 


- 


D 


21.52 


20.84 


— 




.68 


MGPS2 


J082437— 463505 


08 


24 


38 


-46 


35 


09 


263. 


6 


-5. 


1 


2 


3.3 x 10~ 5 


44.3 


33 


.6 


29.6 





.3 





20 


- 


N 


— 


— 


— 




— 


MGPS2 


J083322-444141 


08 


33 


22 


-44 


41 


39 


263 





-2. 


8 


3 


2.0 X 10~ 8 


102.4 


83 


.3 


64.6 





.2 





23 


- 


C 


— 


21.03 


— 




— 


MGPS2 


J083343— 431958 


08 


33 


43 


-43 


19 


58 


261 


9 


-1. 


9 


2 


3.8 X 10~ 6 


26.2 


18 


3 


12.4 





.5 





.36 


- 


N 


— 


— 


— 




— 


MGPS2 


J083420-385401 


08 


34 


20 


-38 


53 


57 


258 


4 


+0 


.8 


3 


3.8 X 10 -6 


66.5 


48 


.8 


40.8 





.3 





24 


33.7 


N 


— 


— 


— 




— 


MGPS2 


J083524— 482342 


08 


35 


24 


-48 


23 


42 


266 


1 


-4. 


7 


3 


4.6 X 10 -7 


51.8 


36 


.9 


31.4 





.3 





.24 


- 


C 


22.12 


— 


19.50 




— 


MGPS2 


J083555-480638 


08 


35 


54 


-48 


06 


38 


266. 





-4. 


5 


3 


4.6 X 10 -6 


41.6 


31 


1 


26.9 





.2 





.21 


- 


c 


19.94 


18.88 


18.31 


1 


.06 


MGPS2 


J090857-523612 


09 


08 


57 


-52 


36 


13 


272 


8 


-3. 


3 


3 


5.7 x 10~ 7 


41.7 


31 


.1 


25.6 





.3 





.24 


- 


c 


— 


— 


— 




— 


MGPS2 


J091052-473743 


09 


10 


52 


-47 


37 


41 


269 


4 


+0 


3 


2 


< 1 x 10~ 14 


34.3 


15 


.6 


10.8 





.7 





.52 


- 


N 


— 


— 


— 




— 


MGPS2 


J095141— 550630 


09 


51 


41 


-55 


06 


31 


279 


1 


-0. 


8 


2 


3.0 X 10 -5 


42.4 


33 


.3 


29.4 





3 





.18 


- 


N 


— 


— 


— 




— 


SUMSS 


J100222— 423135 


10 


02 


22 


-42 


31 


33 


272 


7 


+ 10. 


2 


2 


2.1 x 10 -5 


25.0 


16 


.7 


13.6 





4 





.30 


- 


N 


— 


— 


— 




— 


SUMSS 


J101048— 805633 


10 


10 


47 


-80 


56 


32 


296. 


7 


-20. 


1 


2 


3.7 x 10~ 5 


19.5 


13 


.4 


11.0 





.4 





.28 


- 


N 


— 


— 


— 




— 


MGPS2 


J104530-505429 


10 


45 


30 


-50 


54 


30 


283. 


5 


+7 


.2 


2 


9.6 x 10~ 10 


94.7 


68 


.2 


57.3 





.3 





.25 


- 


C 


— 


— 


— 




— 


SUMSS 


J113907— 462506 


11 


39 


06 


-46 


25 


06 


290. 


2 


+ 14. 


7 


2 


< 1 x 10 -14 


35.7 


18 


.3 


14.3 





.6 





.43 


- 


c 


21.34 


19.97 


19.48 


1 


.37 


SUMSS 


J120525-470714 


12 


05 


26 


-47 


07 


10 


294 


9 


+ 15. 





3 


1.6 X 10~ 5 


29.0 


22 


1 


16.4 





.3 





.28 


- 


N 


— 


— 


— 




— 


MGPS2 


J161340-640022 


16 


13 


40 


-64 


00 


22 


323 





-9. 


3 


2 


9.3 X 10" 6 


34.9 


26 


.9 


23.0 





.3 





21 


- 


C 


20.09 


18.70 


18.36 


1 


.39 


MGPS2 


J170932-441750 


17 


09 


32 


-44 


17 


49 


343 


2 


-2. 


5 


3 


2.1 x 10~ 6 


36.1 


30 


.0 


21.4 





.3 





.26 


- 


c 


20.48 


19.07 


18.52 


1 


.41 


MGPS2 


J171300-541538 


17 


12 


60 


-54 


15 


39 


335. 


4 


-8. 


8 


2 


1.7 X 10~ 6 


33.4 


26 


.9 


20.0 





.4 


o 


.25 




c 


20.76 


20.16 






.60 


SUMSS 


J173100-533721 


17 


31 


00 


-53 


37 


21 


337 


5 


-10 


7 


2 


3.2 x 10 -5 


55.8 


46 


.0 


39.3 





.2 





.17 




c 


20.85 


19.73 


19.13 


1 


.12 


SUMSS 


J175605-464450 


17 


56 


06 


-46 


44 


49 


345 


6 


-10 


7 


2 


1.1 x 10~ 8 


34.1 


22 


.3 


16.8 





.5 





.34 




c 


16.71 


15.56 


15.00 


1 


.14 


SUMSS 


J183056-464744 


18 


30 


56 


-46 


47 


43 


348 


2 


-16. 


1 


3 


5.2 x 10 -6 


16.4 


12 


1 


8.3 





.3 





33 




c 












SUMSS 


J190344-505108 


19 


03 


44 


-50 


51 


08 


346 





-22 


6 


3 


1.5 x 10 -10 


44.9 


33 


.3 


23.1 





.3 





.32 




c 


19.73 


19.25 


19.48 




.48 


SUMSS 


J200936-554236 


20 


09 


35 


-55 


42 


31 


342 


2 


-33 





4 


4.2 x 10~ 6 


21.2 


16 


.8 


12.2 





.3 





.27 




G 












SUMSS 


J201524-395949 


20 


15 


24 


-39 


59 


48 





9 


-32 


5 


2 


2.0 X 10~ 5 


31.5 


23 


.3 


16.4 





.4 





.31 


8.5 


G 












SUMSS 


J202536-813918 


20 


25 


36 


-81 


39 


19 


311 


8 


-30 


2 


2 


1.3 x 10~ 7 


62.1 


47 


.8 


39.7 





.3 





.22 




D 


20.26 


19.88 






.39 


SUMSS 


J211009-772600 


21 


10 


09 


-77 


25 


58 


315 


4 


-33 


7 


2 


1.3 x 10~ 10 


39.0 


25 


.9 


20.6 





.4 





.31 




N 












SUMSS 


J214235-724745 


21 


42 


35 


-72 


47 


40 


318 


9 


-37 


9 


3 


1.8 x 10~ 6 


14.3 


9 


8 


7.8 





.3 





.30 




N 












SUMSS 


J223017-614701 


22 


30 


16 


-61 


47 


03 


327 


1 


-48 


2 


50 


3.0 X 10~ 5 


58.4 


47 


.0 


37.9 





.1 





.21 




D 


21.14 


19.16 


18.02 


1 


.98 


SUMSS 


J223225-615308 


22 


32 


24 


-61 


53 


13 


326 


8 


-48. 


3 


33 


< 1 X 10 -14 


257.7 


221 


.1 


167.5 





.1 





21 




N 












SUMSS 


J225444-500528 


22 


54 


44 


-50 


05 


30 


339 


1 


-58 





2 


3.5 X 10 -7 


24.5 


17 


.0 


12.7 





.4 





31 




E 


20.92 


20.46 






.46 



Table 6. 15 radio sources classified as transient. The columns are: the source name, right ascension and declination (J2000), Galactic longitude and latitude, number of detections, 
number of accepted observations, probability of xf °f th e residuals to a line of constant flux density, maximum flux density at 843 MHz in mjy beam -1 , the maximum and minimum 
detection significance, NVSS flux density at 1.4 GHz in mjy beam" 1 (if covered by NVSS area, ND if not detected), SuperCOSMOS classification: D = point source detected, N = 
no source detected, C = more than 1 point source in the error box, E = extended according to SuperCOSMOS catalogue, G = clearly resolved galaxy, the SuperCOSMOS optical 
magnitudes in blue (B), red (R) and infra-red (I), and the B— R colour. More information on these sources, including identifications, can be found in Table[7] 



Source name 


(X 

(h 


m 


s) 


6 

(d m s 


) 


(de f 


i 
5) 


6 

(deg) 






1 - p(xf c ) 


Smax 
(mjy beam -1 ) 


d ma 




<'i,,i 


n 


NVSS 
(mjy beam -1 ) 


Supcrcos 


B 

(mag) 


R 

(mag) 


I 

(mag) 


B-R 

(mag) 


SN1987A 


05 


35 


28 


-69 


16 


07 


279 


.7 


-31 


.9 


2 


13 


6.5 X 10 -13 


151.8 


13. 


5 


-0 


.4 


- 


c 


18.16 


13.25 


12.65 


4.91 


SUMSS J055712-381106 


05 


57 


12 


-38 


11 


05 


244 


.3 


-26 


.5 


1 


2 


2.5 X 10 -3 


16.7 


7. 


1 


2 


.9 


5.2 


E 


17.89 


16.80 


16.18 


1.10 


J060938— 333508 


06 


09 


38 


-33 


35 


08 


240 


.3 


-22 


.7 


1 


3 


1.4 X 10~ 6 


21.3 





6 


2 


.3 


ND 


G 










J061051-342404 


06 


1(1 


50 


-34 


24 


04 


241 


.2 


-22 


.8 


2 


4 


3.4 x 10 -3 


14.3 


8. 


3 


2 


.2 


ND 


N 










SUMSS J062636-425807 


06 


26 


36 


-42 


58 


07 


251 





-22 


.4 


1 


2 


7.0 X 10~ 4 


19.3 


15. 





2 


.9 




D 


19.39 


18.18 


16.81 


1.20 


J062716-371736 


06 


27 


16 


-37 


17 


36 


245 


.3 


-20 


.5 


2 


3 


3.1 X 10~ 3 


14.9 


6. 


6 


2 


.9 


2.7 


N 










J064149-371706 


06 


41 


49 


-37 


17 


06 


246 


.3 


-17 


.8 


1 


2 


4.4 x 10~ 4 


14.4 


6. 


1 


2 


.4 


ND 


N 










SUMSS J102641-333615 


1(1 


26 


40 


-33 


36 


14 


271 


.3 


+20 


.2 


1 


2 


9.8 X 10 -3 


17.2 


7 


3 


2 


.7 


ND 


N 










SUMSS J112610-330216 


11 


26 


09 


-33 


(12 


13 


283 





+ 26 


.5 


1 


2 


4.6 x 10~ 3 


16.6 


7. 


8 


2 


.8 


7.1 


N 










Nova Muscae 1991 


11 


26 


27 


-68 


40 


31 


295 


.3 


-7 


.1 


1 


2 


< 1 X 10~ 14 


140.4 


67. 


8 





.1 




C 


21.57 




18.72 




J121032-381439 


12 


10 


32 


-38 


14 


39 


204 


.2 


+ 23 


.9 


1 


2 


7.4 x 10~ 3 


15.0 


9. 


.5 


2 


.8 


9.6 


N 










J135304-363726 


13 


53 


04 


-36 


37 


26 


316 


.5 


+24 


.6 


1 


3 


1.0 X 10~ 3 


16.3 


6. 


6 


2 


.9 


5.1 


D 




20.71 






J153613-332915 


15 


36 


13 


-33 


29 


15 


338 


.2 


+ 17 


.9 


1 


2 


1.5 x 10 -3 


18.4 


6. 


2 


2 


.3 


7.7 


N 










GRO 1655-40 


16 


54 


00 


-39 


50 


44 


345 


.0 


+2 


.5 


2 


9 


< 1 X 10~ 14 


6504.2 


17. 


4 


-0 


.8 


ND 


E 


18.58 


15.26 


13.96 


3.32 


SUMSS J224152-300823 


22 


41 


52 


-30 


08 


23 


19 


.1 


-61 


.5 


1 


3 


9.6 x 10 -3 


15.6 


9. 


4 





.6 


ND 


N 
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4.3 Transient sources 



Table [6] lists 15 sources we classify as transients by the cri- 
teria described in £13.61 Eight sources that satisfied these 
criteria are not shown in the table, because visual inspec- 
tion showed incorrect classification due to image artefacts, 
such as nulls causing non-detections or false detections from 
grating rings. 



4.4 Optical Counterparts 

The deepest optical survey of the area of s ky covered by 
our s urvey is the SuperCOSMOS Sky Survey l|Hamblv et al.l 
1200 j ). which was produced by digitising photographic plates 
from the UK Schmidt Telescope. The SuperCOSMOS Sky 
Survey has image and catalogue data in three colours with 
limits down to B ~ 22, R ~ 20 and / ~ 19. To determine 
an optical match for each radio source, we visually inspected 
overlays of the MOST data sources with the SuperCOSMOS 
B image, and searched the SuperCOSMOS catalogues. 

From visual inspection of the overlay plots we were im- 
mediately able to identify four of the 68 transient and vari- 
able radio sources have radio positions directly on, or within 
10 arcsec of the centre of a resolved optical galaxy. These 
sources that have an optical classification of 'G' in Table [S] 
and Table [6] To determine the probability of random false 
association, we visually inspected images of 68 random po- 
sitions on the sky and found that no positions fell near re- 
solved optical galaxies of similar size. We conclude that the 
radio sources that are positionally coincident with resolved 
optical galaxies are very likely to be physically associated. 

Many of our radio sources are found in optically 
crowded regions in the Magellanic Clouds. T he positional 
uncer tainty for MOST sources is < 5 arcsec |Mauch et al.l 
12003! ) ; therefore, we consider any region with 2 or more op- 
tical matches within 5 arcsec to be crowded. Such sources 
have an optical classification of 'C in Table [5] and Table [6l 
and optical magnitudes shown in the tables are those of the 
nearest candidate to the known radio position. 

We consider any radio source with a single optical coun- 
terpart within 5 arcsec to be a detection. These sources have 
an optical classification of 'D' if the match is classified as 
point-like (mean class is 1) in the SuperCOSMOS catalogue, 
or 'E' if the match is extended in the SuperCOSMOS cata- 
logue (mean class is 2). To measure the probability of a false 
match we compare the number of optical sources contained 
within a circle centred on the radio source position, and on a 
random position in the sky (Figure [T7|) over a range of radii. 
We find that, within 5 arcsec, there are 9 optical matches for 
the random positions, and 27 optical matches for the radio 
source positions, corresponding to a false match probability 
at the 5 arcsec match radius of 33 per cent. 

The remaining radio sources have no point-like optical 
counterparts within 5 arcsec to within B ~ 22, R ~ 20 and 
I ~ 19. These sources have an optical classification of 'N' in 
Tableland Tabled 



4.5 Gamma Ray Burst Counterparts 

Some of our transient and variable sources could be related 
to GRB afterglows. In the standard picture, radio afterglow 
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Figure 17. Histogram of the number of SuperCOSMOS B 
sources as a function of radius for randomly chosen fields, and 
for uncrowded MOST transient and variable sources. 



at 843 MHz would be detected several weeks after a GRB 
event and would fade over several months. 

We searched the master catalogue of the Swift GRB 
satellite that has been operational since 2005 Jan 0, for 
GRBs within 10 arcmin of all transient and variable radio 
sources. We found no GRBs within the position uncertain- 
ties of any transient or va riable radio sources . 

We also searched the IStern et alj l|200ll ) catalogue for 
GRBs detected by the BATSE instrument aboard the Comp- 
ton Gamma-Ray Observatory, which was operational be- 
tween 1991 and 2000, for our transient sources only. For each 
transient source, between 5 and 10 counterparts were identi- 
fied within positional uncertainties (which are typically tens 
of degrees) . No GRBs were detected in the weeks before the 
first detection of a transient radio source with a typical GRB 
afterglow light curve. 

The vast majority of GRBs are detected at z > 0.1 
|jakobsson et all l2006h with the remainder being detected 
in gamma rays at much larger distances. All of our transient 
sources with resolved optical counterparts are at z < 0.1. 
Hence, if our transient sources had a gamma-ray trigger, we 
would expect them to be detected if a GRB satellite had 
been operational at the time. 

We note that a number of transient sources with 
radio light curves consistent with GRB afterglows (i.e. 
J060938-333508, J062636-425807, J064149-371706 and 
SUMSS J102641-333615), would not have had a GRB de- 
tection in either Swift BAT or BATSE catalogues as neither 
instrument was fully operational when the first MOST de- 
tection was made. 



5 DISCUSSION 
5.1 Overall statistics 

Our sample comprises light curves for 29730 sources, with 
a two-epoch equivalent sky coverage (c.f. Equation [1]) of 

3 http : //www. swift . ac.uk/swift_live/arnie5.php 
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Figure 18. Log N-Log S plot for the combined SUMSS and 
MGPS-2 catalogue, all sources used for our variability analysis, 
and the sources classified as variable and transient. The error bars 
are Poisson errors. The sources available for variability analysis 
are incomplete above 100 mjy beam -1 , as strong sources have 
poor fits in MOST images due to artefacts. 



2775.7 deg . The shortest interval between accepted mea- 
surements of a single source is 24 h and the longest interval 
between accepted measurements is 7301 days (~ 20 years). 

For a source to be in our sample, it must appear in 
at least two images for which the exact two-epoch coverage 
applies, rather than the equivalent coverage. The two-epoch 
(exact) coverage is approximately 2000 deg 2 (c.f. Figure [3} 
resulting in a source density of 15 deg" 2 . The SUMSS source 
density is 31.6 deg -2 meaning the overall completeness of 
our sample is 47 per cent. 

A Log N-Log S plot for the SUMSS and MGPS-2 cat- 
alogues, our whole sample, the transient sources and the 
variable sources is shown in Fig. [18] At low flux densities 
the catalogue is incomplete due to variations in the noise 
floor and at high flux densities it is due to artefacts around 
bright sources. In general the low completeness is due to the 
strict requirements on goodness of fit we imposed in order 
to remove false positives due to artefacts. 

Fig. \W\ and [2S1 show histograms of modulation index, 
and V for all sources in our sample. 



5.2 Variability time-scales 

If there is a sufficient number of measurements on the 
light curve of a source, a structure function can be formed 
to make an estimate of t he variability tim e-scale (e.g. 
iGaensler fc HunsteacfeOOOl or lLovell et "atboOSl ) . As the ma- 
jority of our sources have only a small number of measure- 
ments (see Fig. Hip , structure functions cannot be used. To 
determine if there is a preferred variability time-scale for 
our sources, we compare the distribution of intervals be- 
tween measurements of our variable sources, with that of 
the whole sample as follows: For each source, we take each 
non-redundant pair of flux density measurements and bin 
by the time interval between these measurements, i.e. if 
there are iV mea s measurements on the light curve, there are 
Nmeas (AVeas — l)/2 intervals. We then form a histogram 
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Figure 19. Histogram of modulation indices (m = cr/S, see 
£|3. 5.21 1 for all sources, including variable sources. 
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Figure 20. Histogram of fractional variability (see Equation [8J 
for all sources and variable sources. 



comprised of all sources, and the subset of variable sources 
in Table[5]only. If the variable sources have a preferred time- 
scale, this should be evident in the variable sources consti- 
tuting a larger fraction of their measurement intervals at a 
particular time-scale. Fig. [5] depicts the histograms of inter- 
vals described above. 

Fig. [21] shows the ratio of the number intervals belong- 
ing to the variable and transient sources to the number of 
intervals belonging to the whole sample peaks. The ratio 
peaks at around 1000 days and then falls rapidly by 2000 
days. This suggests that any blind survey for variable radio 
sources will detect the majority of variable sources in the 
first 1000 to 2000 days. 

5.3 Source counts for variable sources 

The number of variable sources in our sample can be deter- 
mined in several ways. The source list in Table [5] contains 
55 sources and we are confident that at least 90 per-cent of 
them are truly variable, which implies a source count of 0.020 
sources per two-epoch square degree, but this is not the 
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Figure 21. The ratio of the number of non-redundant pairs of 
intervals in the light curves for for variable and transient sources, 
to the total number of intervals for all sources, as a function of 
the interval in deays. The y errors are the Poisson errors. 



whole story. Fig. [13] shows there is a clear excess of sources 
with P(xi 2 c ) < 0.1. While we cannot determine whether a 
particular source's variation is attributable to measurement 
errors or true variation, we can compare the number of ap- 
parently variable sources to the expected number if only 
Gaussian measurement errors were present. If we conserva- 
tively take the 'expected' line as the number of non-varying 
sources, we arrive at an excess of 200 sources, corresponding 
to a lower limit on the the two-epoch areal density of vari- 
able sources at 843 MHz of 0.072 deg -2 . This would appear 
to be an underestimate, as the actual number of sources with 
P(Xio) < 0.2 is substantially less than the expected number, 
due to over-estimated errors. If we extrapolate the actual 
number of sources from P(xic) ~ 0-2, we arrive an an excess 
of approximately 350 sources with -P(xfc) < 0.2 at 843 MHz 
with flux densities between 15 and 100 mjy beam -1 , which 
corresponds to a two-epoch areal density of variable sources 
at 843 MHz of 0.126 deg -2 . After accounting for complete- 
ness of 47 per cent, the areal density rises to 0.268 deg -2 in 
2 epochs. 

A definitive statement about the number of variable 
sources requires a caveat on the time-scale of the variability. 
Our data are not regularly sampled and so we can only make 
broad statements about variability on time-scales greater 
than 1 day and less than about 10 years. 



5.4 Interpretation as interstellar scintillation 

At 843 MHz, our observations are firmly in the strong or 
refra ctive reg i me of interstellar scintillation. In these condi- 
tions iRickettJ (1 19861 ) derives the maximum rate of change of 
flux of sources with a brightness of T\i x 10 12 K at latitudes 
|6| > 10° as 



[dS/dt) m ax < Srms/r ~ 30Ti 2 (sin 6) ' 5 mjy day 1 



(10) 



giving a range of (dS/dt) ma x for a source with a brightness 
temperature at the inverse Compton limit of between 12.5 
and 30 mjy day -1 . As our data are sparsely sampled in 



Figure 22. Fraction of variable and transient sources vs Galactic 
latitude. 



the time domain, this limit is easily satisfied for all of our 
sources. 

The number of sources that scintillate is expected to 
increase as the line of sight to the source becomes closer to 
the Galactic plane, and is more likely to pass through a scin- 
tillating screen of material. In Figure[22]we plot the fraction 
of sources that we defined as variable and transient, binned 
by Galactic latitude. The plot suggests, at best, a mild de- 
pendence on Galactic latitude. Nonetheless, we expect that 
refractive interstellar scintillation of compact radio sources 
is a likely explanation for a large fraction of our v ariable 
sources, as also seen bv lGaensler fc Hunsteadl (I2000T ). 

Some 30 per cent of the variable sources have no op- 
tical counterpart. In the absence of more data, we cannot 
positively identify these sources. 



5.5 Source counts for transient sources 

Transient emission on characteristic time-scales of less than 
the full synthesis time is unlikely to be detected by our tech- 
nique, because a significant change in flux density would re- 
sult in a distorted synthesised beam for that source, which 
would result in large errors in a Gaussian fit and ultimately 
rejection from the analysis. Therefore, we will omit from our 
discus sion source classes su ch as giant pulses from neutron 
stars (Hankins ct al. 2003) as well as flares from extraso- 
lar pla nets dBastian et al„ 2000) and fast extragalactic tran- 
sients (|Lorimer et al.l 2007h . which have time-scales much 
shorter than 12 h. Also, a source that is constant over 12 h 
but bright (> 100 mjy), or in a region containing bright 
sources (such as the Galactic plane), is unlikely to be de- 
tected due to image artefacts or large formal errors in source 
fitting. 

Table [6] shows 15 sources that satisfied our definition 
of a transient with a flux density above 14mJy beam - , 
corresponding to an areal density of transient sources at 
843 MHz of 6 x 10 -3 deg -2 in 2 epochs. After accounting 
for completeness of 47 per cent, the areal density rises to 
1.3 x 10 -2 deg -2 in 2 e pochs. 

iBower et all ll2007h compute a two-epoch areal density 
of transient sources of 1.5 deg -2 above 370 /iJy at 5 GHz 
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and a flux density distribution following S^ 1 with 7 = 1.5, 
which corresponds to an areal density of 2.0 x 10~ 2 deg" 2 , 
assuming a spectral index of -0.7 and flux density limit of 
15mJy beam -1 , in broad agreement with our findings. We 
note that there are aspects of our transients which are not 
the same as those described bv lBower etTafl (|2007l ). In par- 
ticular, our transients are detected at a frequency an order of 
magnitude lower, have flux densities an order of magnitude 
higher, and have upper limits on characteristic time-scales 
of one or two orders of magnitude longer. 

5.5.1 Known transients 

Our search detected three previously identified radio tran- 
sients: 

(i) GRO 1655-40 (J165400-395044) is a low mass X-ray 
binary (or microquasar) , detected in our survey because of 
the ima ges taken while monitoring a drama tic radio flare in 
1994 (|Campbell- Wilson fc Hunsteadlll994ah . 

(ii) Nova Muscae 1991 (J112627-684031) is a low mass 
X-ray binary, detected in ou r survey from an image taken 
during the outburst in 1991 jKesteven fc Turtle|[T99ll 'l. 

(iii) SN1987A (J053528-691607) was a supernova in the 
Large Magellanic Cloud. It is det ected in our sur vey due 
to a regular monitoring campaign |Ball et al.ll200lf ) during 
which its flux density has been steadily increasing. 



6 DISCUSSION OF SELECTED SOURCES 

In this section we discuss sources with noteworthy light 
curves or optical counterparts. 

6.1 Variable sources with resolved optical 
counterparts 

6.1.1 SUMSS J011019-455112 

The MOST light curve of SUMSS J011019-455112 fFig. 123) 
shows two decreasing measurements over 10 years, followed 
two years later by four significantly higher measurements 
closely spaced in time. We have also obtained flux density 
and position measurements from archival ATCA observa- 
tions at 1.4 GHz used for the Phoenix Deep Field survey, 
which rep orts the average flux density over all epochs of 
10.6 mJy (|Hopkins et alj|2003h . The optical overlay of the 
higher resolution ATCA data shows the source is centred 
9.6 arcsec from the centre of the extended 2MASS source 
2MASX J01101993-455118, an SO galaxy in the gala xy clus- 
ter Abell 2877 at z = 0.023 (|Caldwell fc Roselll997l ). Using 
the highest MOST flux density, the implied isotropic radio 
luminosity is L v — 2 x 10 29 erg s" 1 Hz -1 . The spectral index 
between the mean 1.4 GHz flux density and the minimum 
and maximum flux densities from MOST is in the range to 
-1.3 but we note that no measurements are contemporane- 
ous. The radio source detected by ATCA at 1.4 GHz was a 
point source at the 12 arcsec resolution of the observations, 
and had a position consistent with a = lhlml9.442s S = 
-45d 51ml3.958s (J2000) in all ATCA epochs. 

As the MOST contours are offset from the optical centre 
of 2MASX J01101993-455118, it is unlikely to be associated 



with AGN activity, so a Radio Supernova (RSN) or GRB 
afterglow are possible candidates. 

Radio Supernovae (RSNe) have typical radio spectral 
luminosity of L v ~ 10 25 erg s" 1 Hz" 1 , but can re ach as high 
as L v ~ 2 x 10 27 erg s" 1 Hz" 1 (|Eck et al.l [20021 ') and have 
fading times of 20 years or more. SUMSS J011019-455112 
has a much higher radio luminosity than even the most lu- 
minous known RSN (1979C), so if it is a radio supernova, it 
is unusually bright. 

The radio luminosity of SUMSS J011019-455112 is 
closer to that of GRB afterglows, which typ ically have a 
range of L v ~ 6 x 10 28 - 4 x 10 31 erg s" 1 Hz" 1 (|Weiler et all 
2002). On the other hand, the 12 year time-scale is typi- 
cal for a Type II RSN but too long for a GRB afterglow, 
which typically fade after several weeks. The approximately 
constant flux density for the MOST and ATCA observations 
could be explained by pre-existing disc synchrotron emission 
from the host galaxy, with the flare in the MOST observa- 
tions being due a flare or stellar explosion overlaid. In such 
a case, one would expect the disc synchrotron emission to 
appear resolved at the ATCA resolution and roughly fol- 
low the galaxy optical emission, but we find no evidence for 
resolved emission in the ATCA observations. 

We conclude that SUMSS J011019-455112 is associ- 
ated with an unusual stellar event or explosion in 2MASX 
J01101993-455118. 



6.1.2 S UMSS J 201 524 - 39594 9 

The light curve for SUMSS J201524-395949 (Fig. 121) shows 
an NVSS detection and two MOST detections over 10 years. 
The MOST flux density doubled over 4 years. The MOST 
contours are centred on a barred spira l galaxy ESO 340- G06 
at redshift z = 0.02 (|jones et al.ll2004h implying, if the radio 
source is associated with the spiral galaxy, an isotropic radio 
luminosity of L y = 3 x 10 29 erg s" 1 Hz" 1 . 

The 6dFGS spectrum l| Jones et al.l 12004^ 1 shows strong 
Balmer emission lines, indicative of a star-forming galaxy. 
The IRAS 60 /im fl ux is 1 Jy, but the galaxy is not detected 
at 12 and 25 /xm (|Moshij[l990l 'l. suggesting the presence 
of cold dust. It is rare, but not impossible, for such a star- 
forming galaxy to contain an AGN. 

The radio flux in this case is likely due to a stellar ex- 
plosion, as widespread star formation or disc synchrotron 
emission cannot explain the radio variability. Once again, 
the time-scale of 10 years argues strongly against a GRB af- 
terglow. RSN light curves are known to behave differently at 
different frequencies, with flux density increase being later 
and slower at lower frequencies. The increasing flux density 
over 4 years at 843 MHz is consistent with a Type II RSN 
but the 6 year interval between MOST and NVSS detec- 
tions is arguably too long to be due entirely to an RSN, in 
which case it may be that the RSN occurred after 1996 and 
the NVSS and first MOST measurement could be attributed 
to pre-existing disc synchrotron emission or widespread star 
formation. 

6.1.3 S UMSS J 200936- 554236 

The light curve for SUMSS J200936- 554236 (Fig. [25} shows 
an increase in flux density over 3 years, and then steady 
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Figure 23. Variable source SUMSS JOlf Of 9 -455U2. Top panel : 
Radio contours from the Phoenix deep field (Hopkins ct al. 2003) 
at f.4 GHz at 2 to 8 mjy beam -1 in steps of 2mjy beam -1 
overlaid on a SupcrCOSMOS B image. Bottom panel: Radio light 
curves at f.4 GHz (ATCA archive) and 843 MHz (MOST). 
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Figure 24. Variable J20f 524-395949. Top panel: MOST con- 
tours from 5 to 30mjy beam -1 in steps of 5mjy beam - 1 for the 
brightest epoch (2004 May Of) overlaid on a SuperCOSMOS B 
image. Bottom panel: light curve showing an NVSS detection and 
a doubling in 843 MHz flux density over 5 years after a 6 year 
interval. 



flux density over 2 years. The MOST contours are centred 
on the barred spira l galaxy IC4957 at redshift z = 0.032 
l|Strauss et all 1 19921 ). This galaxy is detected across the 
IRAS bands at 0.2 Jy (25 fim), 2.4 Jy (60 fim) and 4.5 Jy 
(100 /irn), implying warm dust. The implied radio luminos- 
ity, if the radio source is associated with the spiral galaxy, 
is L u = 5 x 10 29 erg s" 1 Hz" 1 . 

As the MOST contours are centred on the optical 
source, AGN activity or ISS of an AGN is a possible 
explanation but AGNs in barred spiral galaxies are rare 
ijWilson fc Colbert! 1 19951 ). The warm dust implies star for- 
mation, so a stellar event is conceivable. The 3 year time- 
scale argues strongly against a GRB afterglow interpre- 
tation, and the shape and time-scale are consistent with 
a Type II RSN. The luminosity, once again, is higher 
than that of any known RSN. We conclude that SUMSS 
J200936— 554236 is associated with an unusual stellar event 
or explosion in IC4957. 



6.2 A compact steep spectrum source: SUMSS 
J223225-615308 

The light curve for SUMSS J223225-615308 (Fig. [26) shows 
two pronounced dips in flux density over two years and it 
is the most dramatic of some five variable sources in this 
field. There are other sources at the same flux density level 
in this field that do not exhibit variability, ruling out a cal- 
ibration effect. An overlay of SUMSS J223225 -615308 and 
the P arkes-MIT-NRAO (PMN) survey data (|Wright et all 
1 19941) confirms an association with PMNJ2232— 6153, whose 
reported flux density is 40 ± 7mJy at 4.85 GHz. This flux 
density together with the range of flux densities shown in the 
843 MHz light curve, implies a spectral index between —0.8 
and -1.0, although we note that the PMN and MOST ob- 
servations are not contemporaneous. If the variation is due 
to scintillation, this implies a compact source or component, 
and together with the spectral index implies a classification 
as a compact steep spectrum source. 
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Figure 25. Variable SUMSS J200936-554236. Top panel: MOST 
contours at — 1.5mjy beam -1 and 6 to 18mjy beam -1 in steps 
of 3 mjy beam -1 for the brightest epoch (1999 March 15) overlaid 
on a SuperCOSMOS B image. Bottom panel: radio light curve 
indicating a rising and then constant flux density. 



6.3 Transient sources with resolved optical 
counterparts 

6.3.1 J 060938- 333508 

The light curve for J060938-333508 (Fig. [27]) shows MOST 
and NVSS non-detections followed by a single MOST de- 
tection. Unlike the vast majority of sources in the MOST 
archive, the MOST contours appear rotated with respect to 
the MOST beam possibly indicating a change in flux den- 
sity over the 12 hr synthesis time. The MOST contours are 
centred 9 arcsec from the centre of Fairall 1138. a g alaxy 
with spectral type SBab D dDressler fc Shectmarj|l988h and 
with redshift z = 0.037 ([Dale et allll998T ). 

Assuming the radio source and galaxy are associated, 
the inferred isotropic radio luminosity from the brightest 
epoch (2004 December 9) is L v ~ 6 x 10 29 erg s" 1 Hz -1 at 
843 MHz. 

The offset from the centre of the optical galaxy, and 
the fact that spiral galaxies rarely contain an AGN argues 
against an AGN source for the radio variability. The optical 
morphology of Fairall 1138 appears quite distu rbed and it 
has a 60 /im flux of 345 mjy (Mo shir et alj|l992h . which im- 
plies the presence of cold dust. The disturbed morphology, 
spectral type and presence of dust in Fairall 1138 support 
the interpretation that it is undergoing star formation, so a 
stellar explosion, either GRB or RSN are possible explana- 
tions. 

The spectral luminosity is very high for a RSN and we 
are unable to discriminate between Type Ib/c or Type II 
RSNe by the light curve time-scales, as the time interval 
between the detection and non-detection epochs is too large. 

The spectral luminosity of J060938— 333508 is within 
the range of GRB afterglows. It is also within the error cir- 
cle of GRB 940526B, which occurred after the MOST non- 
detection in 1993, and 10 years before the MOST detection. 
If J060938-333508 is the radio afterglow of GRB 940526B, 
then the radio detection 10 years after the gamma ray event 
is unlike known GRB afterglows, which typically peak at 
843 MHz a few weeks after the explosion and fade over about 
3 years. We consider an association of GRB 940526B and 
J060938-333508 unlikely. 

We consider an unusual stellar event in Fairall 1138 as 
the likely interpretations of this source. 
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Figure 26. Light curve of the candidate compact steep spectrum 
source: SUMSS J223255-615308 



6.3.2 SUMSS J055712-381106 

The light curve for SUMSS J055712-381106 (Fig. [28} shows 
an NVSS detection followed by a MOST detection approx- 
imately 10 years later and then a non-detection 6 days af- 
ter that, consistent with either a flaring source or a highly 
variable source occasionally appearing above our sensitivity 
limit. The MOST contours appear slightly elongated and ro- 
tated with respect to the MOST beam possibly indicating 
a change in flux density over the 12 hr synthesis time. Flar- 
ing or scintillating AGN or flaring radio stars are possible 
counterparts with these properties. 

The SuperCOSMOS B image shows what appears to be 
a blend of three objects, a star-like object to the south, a 
faint star-like object immediately to its north, and an ex- 
tended object to the north-east. The MOST and NVSS ra- 
dio sources cannot be conclusively associated with any of the 
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Figure 27. Radio transient J060938-333508. Top panel: MOST 
contours from 6 to 18mjy beam -1 in steps of 3mjy beam -1 
for the brightest epoch (2004 December 9) overlaid on a Supcr- 
COSMOS B image. The MOST synthesised beam is drawn at 
the bottom left. Bottom panel: MOST and NVSS light curves 
indicating a single detection in 2004 at 21.3 mjy beam -1 and 3 
non-detections . 



Figure 28. Radio transient SUMSS J055712-381106. Top 
panel: MOST contours from 3 to 15mjybeam -1 in steps of 
3mjy beam -1 for the brightest epoch (2004 November 7) over- 
laid on a SuperCOSMOS B image. Bottom panel: MOST light 
curve indicating an NVSS and MOST detection, with a non- 
detection 6 days after the MOST detection. 



three sources in the optical blendwhich may in fact be a sin- 
gle interacting system. A more accurate position is required 
for a more definitive association. 

With no clear optical association we cannot positively 
classify this source, but the time-scales between the NVSS 
detection and MOST detection and non-detection rule out 
RSNe and GRB afterglows. 

6.4 Transient sources with point-like optical 
counterparts 

6.4.1 S UMSS J 062636- 4 25807 

The light curve for SUMSS J062636-425807 (see Appendix 
[B} has a detection and a non-detection separated by almost 
a year. The optical counterpart is point like and the colours 
are red, with B-K = 5. The source is most likely an AGN 
scintillating above our noise threshold. 



6.4.2 J 135304- 363726 

The light curve for J135304-363726 (see Appendix [Bj 
is consistent with a flaring source, or a variable source 
occasionally appearing above our sensitivity limit. It is 
marginally detected in the SuperCOSMOS R plates and not 
detected in B or I, or any 2MASS images. This source is 
likely an optically faint AGN scintillating above our noise 
threshold. 



6.4-3 Unidentified variables 

Some 25 per cent of the variable sources have a clear point- 
like or marginally resolved optical counterpart most of which 
have red colours, with < B — R < 2. The optical point 
source and radio detection would typically imply AGN, but 
the colour distribution is unusua l for AGN, for whi ch the 
vast majority have B - R < 1 (|Croom et all Ho04) . The 
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reddest sources are therefore possibly flaring M-dwarfs, but 
optical spectra are needed to confirm this. 

6.4-4 Unidentified transients 

The spectral luminosity of our transients is L v ~ 
10 19 (D_l/1 kpc) 2 erg s _1 Hz _1 , where Dl is the luminosity 
distance. This spectral luminosity i s an order of magn itude 
more luminous than those found bv lBower et~ai1 l|2007t ) at 5 
and 8 GHz, but roughly consistent assuming a spectral in- 
dex of a — —0.7. The sampling interval of our lig ht curves is 
typica lly much longer than those discussed bv lBower et al.l 
(2007) (typically years as opposed to weeks) and we have not 
yet obtained dedicated optical follow-up of our sources, so 
our limits on optical counterparts are not very strong. The 
scarcity of data means we are unable to robustly classify a 
large fraction of our transient sources. 

Of the unidentified transients, none have X-Ray coun- 
terparts in the ROSAT X-ray All-Sky Survey Bright Source 
Catalogue l|Voges et al.lll999l ). GRB afterglows without a 
gamma ray trigger (the so-called orphan GRB afterglows) 
are a possibility and might explain the higher average lu- 
minosity. Radio supernovae are also a possibility as long as 
the hosts are nearby and intrinsically faint, but distant RSN 
would appear to be ruled out by the high luminosity. Stellar 
events are less likely, but flares from nearby late type stars 
are not ruled out. Soft gamma repeaters and X-Ray binaries 
are also possibilities, although the rarity of such objects and 
our measured transient rate are incompatible. 

The most likely source of variability from extragalac- 
tic sources at 843 MHz is refractive Interstellar Scintillation 
(ISS) of the compact components of AGN. By our defini- 
tion, our transients could well be such scintillating sources 
occasionally appearing above our flux threshold, although 
the light curves of the transient sources have V > 0.6, 
which would make them the among most extremely variable 
sources in our sample (c.f. Figure [217)) . Star- forming galax- 
ies, whose radio flux is from widespread star formation and 
diffuse synchrotron emission in the disk, are unlikely to be 
intrinsically variable on our time-scales (due to light travel 
time constraints), or extrinsically variable (due to larger an- 
gular size quenching the interstellar scintillation). 

Another explanation for the transient sources is time- 
integrated emi ssion from p ulsars with on-times of greater 
than ~ 12 h. lOfek et al.l (|2010l ) prop osed pulsars as the 
sour ce of the trans i ents d iscovered by I Bower et al.l (|2007l ) 
and iNiinuma et al.l |2007h based on population arguments, 
and there is some observational evide nce for this. Intermit - 
tent pulsars such as PSR B1931+24 (|Kramer et al.ll2006h , 
have active periods of a few days and then turn completely 
off. PSR J0941— 39, which was discovered as a Rotating 
Radio Transient (RRAT) with 5 single pulses in one ob- 
servation, exhibite d standard pulsar emission during an- 
other observation (|Burke-Spolaor fc Bailcs 2010). Such ob- 
jects might appear as transients in our survey if the time 
between mode-changes had a time-scale of 12 h or longer, 
and the time-integrated flux densities in both modes were 
compatible with our non-detection thresholds. Pulsars are 
also not likely to be optically detected at the SuperCOS- 
MOS plate limits, which would explain the lack of optical 
detections of some our transient sources. 

Finally, microlensing is an unlikely explanation for our 



unknown transients, as the implie d number of sourc es to 
explain our sample is far too high |Bower et al.ll2007r i. Re- 
flected solar flares are unlikely, as a 1 MJy solar flare would 
need to reflect of a 1000 km object at a distance of only 
7000 km to be detected above our threshold. 



Table 7. Summary and classifications of the 15 radio sources classified as transient. 



Source 




Radio Lightcurve 


Counterparts at other wavelengths 


Classification 


&jn 195 ( a 




Linear increase 


Bright point source in SupcrCOSMOS 


Supernova. 








R and I bands. 




oUlvloo JUOOflZ- 


QOI 1 Oft 


Detection followed by a non-detection 


Blend of three objects in SupcrCOS- 


Unknown 






6 days later. 


anu ZMAOO. 




JUDUyoo— ooooU I 




Single detection. Non-zero PA possibly 


9 arcsec from the optical centre of 


t r r\29 1 1 r> -ui u • i~ 

Luminosity L v ~ 6 x 10 erg s Hz . Possible ultra-bright 






indicates variation during 12h synthe- 


Fairall 1138, a barred spiral galaxy at 


KolN or uKd aitcrglow. 






sis. 


Z — U.Uo f . 




Tnci nf:i Q/io/in/i 
JUOlUol — o4z4U4 




Non-detection between two detections. 


No known counterparts. 


TT 

Unknown 


DUlVLoo JUDzOoO— 


A r )KQI~\'7 

-4zooU f 


Single detection. 


oupcroOoMvJo ana zMAoo colours 


Probable scintillating AGN. 








typical of a quasar. 




J062716-371736 




Non-detection between two detections. 


No known counterparts. 


TT 

U nknown 


J064149-371706 




Single detection. 


No known counterparts. 


Unknown 


SUMSS J102641- 


-333615 


Single detection. 


No known counterparts. 


Unknown 


SUMSS J112610- 


-330216 


Detection and non-detection separated 


No known counterparts. 


Less likely to be an AGN or explosion. Possible flare star other 






by 1 day. 




short term variable. 


Nova Muscae 1991 


Single detection. 


Bright point source visible only in ou- 


Low mass X-ray binary. 








perL/ColvlUo ix image, crowded ncld. 










Point source in RASS. 




J121032-381439 




Non-detection between two detections. 


No known counterparts. 


TT 

Unknown 


J135304-363726 




Non-detection between two detections 


T7V. • 4- / ]\/f Of! A Jl/f OH 1\ nnint 

rainc- yivijn — zu.^i — ^u. / y point 


Probable scintillating AGN but faint, red detection could be 








source m oupcrovjoivivjo. 


a star. 


J153613-332915 




Non-detection between two detections. 


No known counterparts. 


Unknown 






Non-zero PA possibly indicates varia- 










tion during 12h synthesis. 






GRO 1655-40 




Power law decrease over 3 months, 


Point source in all SupcrCOSMOS 


GRO J1655-40. A low mass X-ray binary and microquasar. 






then non-detection after 11 years. 


plates. 




SUMSS J224152- 


-300823 


Single detection 


6 arcsec from a blue point-like optical 


Probable scintillating AGN if optical association is correct, 








source. Possible association. 


otherwise unknown. 
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7 CONCLUSIONS 

We have conducted a 22 year survey for radio variability at 
843 MHz and measured 29230 radio light curves. We have 
discovered 15 candidate transient sources and 53 candidate 
highly variable sources. Only 3 of the transients were previ- 
ously known. 

We conclude that many variable sources can be ex- 
plained as scintillating AGN, and some that are associated 
with nearby galaxies may be over-luminous radio supernovae 
with atypical light curves. We conclude that at least 3 of the 
transients are unlike any known source and could belong to 
the class o f radio transients w ithout optical counterpart dis- 
covered bv lBower etafl l|2007h . which have a number of pos- 
sible explanations including giant M-dwarf flares, or flaring 
Galactic neutron stars. 

We have also presented a number of statistical tech- 
niques to aid in future radio variability surveys. In particu- 
lar, we have described techniques to remove systematic gain 
errors, to confirm errors in flux density are correctly esti- 
mated and to choose a robust variability threshold. 

Upcoming wide-field radio variability surveys will offer 
a new insight into this parameter space. We have shown that 
the radio sky is indeed variable so these surveys are likely 
to uncover many new sources. We have also confirmed that 
there is at least one class of transient radio source without 
optical counterparts. The time interval between radio activ- 
ity and optical measurement in our cases was several years, 
highlighting the need for prompt multi-wavelength follow- 
up. 
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APPENDIX A: LIGHT CURVES OF VARIABLE 
SOURCES 

This appendix contains radio light curves of all variable 
sources shown in Table in RA order (Fig.[M|. 843 MHz 
flux density measurements from MOST are shown as black 
squares. Flux density measurements or limits from NVSS, 
where available, are shown as open circles with x er- 
rors indicating the observation duration of the NVSS sur- 
vey. For SUMSS J011019-455112, flux density measure- 
ments at 1.4 GHz are shown as open triangles. For SUMSS 
J223225-615308 the Parke s-MIT-NR AO (PM N) flux den- 
sity measured at 4.85 GHz |Wright et al.|[l994T ) is shown as 
an open triangle. 



APPENDIX B: TRANSIENT SOURCES 

This appendix contains radio light curves of all variable 
sources shown in Table(6] in RA order (Fig. lBl) ). Each source 
is shown on a single row. The leftmost panel is the radio light 
curve. 843 MHz flux density measurements from MOST are 
shown as black squares. Flux density measurements or limits 
from NVSS, where available, are shown as open circles with 
x errors indicating the observation duration of the NVSS 
survey. The centre and right panels show the images of the 
source at the minimum and maximum detection significance, 
in time order from left to right, with the position of the 
source encircled. The grey scale of all images is linear from 
—5 to 20mJy beam - . 

This paper has been typeset from a TgX/ DTpX file prepared 
by the author. 
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Figure Bl. Transient sources. Left panel: radio light curve. The centre and right panels show the images of the source at the minimum 
and maximum detection significance, in time order from left to right, with the position of the source encircled 
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Fig. IBll (continued) . 
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Fig. [Bl] (continued). 
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